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ABSTRACT 



Interest in bat-borne diseases and parasites has grown in the past decade over concerns for human 
health. However, the drivers of parasite diversity among bat host species are understudied as are the links 
between parasite richness and emerging risks. Thus, we aimed at exploring factors that explain macro 
and microparasite species richness in bats from Southeast Asia, a hotspot of emerging infectious diseases. 
First, we identified bat species that need increased sampling effort for pathogen discovery. Our approach 
highlights pathogen investigation disparities among species within the same genus, such as Rhinolophus 
and Pteropus. Secondly, comparative analysis using independent contrasts method allowed the identifi- 
cation of likely factors explaining parasite and viral diversity of bats. Our results showed a key role of 
bat distribution shape, an index of the fragmentation of bat distribution, on parasite diversity, linked 
to a decrease for both viral and endoparasite species richness. We discuss how our study may contribute 
to a better understanding of the link between parasite species richness and emergence. 
© 2014 The Authors. Published by Elsevier Ltd. on behalf of Australian Society for Parasitology. This is an 
open access article under the CC BY-NC-ND license (http://creativecommons.Org/licenses/by-nc-nd/3.0/). 



1. Introduction 

Incidence of emerging infectious diseases (EIDs) has dramati- 
cally increased in recent decades (Jones et al., 2008). A majority of 
EIDs are zoonoses and most of them originate in wildlife (zoonotic 
spillover). EID events are often due to environmental alteration, 
including agricultural intensification and habitat modification, 
global trade and travel (Cunningham, 1996; Daszak et al., 2000). 
Several studies and reviews have suggested that the risk of disease 
transmission from wildlife to humans should increase with biodi- 
versity loss and the expansion of human populations (Dobson 
et al., 2006; Keesing et al., 2009), as humans will get into contact 
with a large pool of known and unknown zoonotic pathogens 
from wildlife (Mahy and Brown, 2000; Murray and Daszak, 2013). 
Bacteria and rickettsia represent the majority of EIDs, viral and 
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prion pathogens cause 25.4% and other parasites 20.3% of EID 
events (Jones et al., 2008). 

Among known reservoir species of viral EIDs, bats, which repre- 
sent 20% of mammal species (Simmons, 2005), play an important 
role in the maintenance and spread of various viral diseases 
(Sulkin and Allen, 1974; Ghatak et al., 2000; McColl et al., 2000; 
Olival et al., 2012; Luis et al., 2013), including members of 
the alphaviruses, flaviviruses, paramyxoviruses, rhabdoviruses, 
coronaviruses and arenaviruses among others. Several notable 
bat-borne viruses are a public health concern in Southeast Asia 
(SEA) including Nipah virus, lyssaviruses and Severe Acute Respira- 
tory Syndrome coronavirus (Mackenzie et al., 2003). While bats are 
increasingly viewed as a threat to human health, these mammals 
have important roles in ecosystems such as pollination, seed dis- 
persal and predation on insects. They are excellent bioindicators 
of environmental changes as they are sensitive to a wide range of 
anthropogenic disturbances such as urbanization, agricultural 
intensification, habitat loss and fragmentation (Clarke et al., 
2005; Jones et al., 2008). Understanding the ecology of these 
potential reservoirs of zoonotic pathogens is needed for improving 



http://dx.doi.org/10.101 6/j.ijppaw.2014.06.003 

2213-2244/© 2014 The Authors. Published by Elsevier Ltd. on behalf of Australian Society for Parasitology. 
This is an open access article under the CC BY-NC-ND license (http://creativecommons.Org/licenses/by-nc-nd/3.0/). 



162 



N. Cay et al. / International journal for Parasitology: Parasites and Wildlife 3 (2014) 161-170 



management of bats and their habitats, ultimately ensuring the 
health of humans, livestock and wildlife species, while keeping 
their functional roles in the ecosystems (Breed et al., 2006). 

The potential of human-bat interactions, either direct or indi- 
rect, may be underappreciated and greater than expected for some 
species. While bats rarely seek direct contact with humans, they 
often roost in or near human dwellings, which can lead to acciden- 
tal contact or exposure to bat excreta. For example, Nipah virus has 
been transmitted to people who have eaten bat-contaminated date 
palm sap (Luby et al., 2006). 

Wild animals are known to host different pathogens at the same 
time (Petney and Andrews, 1998; Drake and Bundy, 2000; Bordes 
et al., 2008; Bordes and Morand, 201 1 ). Polyparasitism (or multiple 
infections) and infection dynamics in host species are important 
features to understand the mechanisms of E1D; the risk of disease 
transmission may depend on the pathogen richness found in natu- 
ral reservoir species (Wolfe et al., 2005). Morse (1993) termed this 
pathogen diversity in wildlife hosts the 'zoonotic pool'. In the liter- 
ature, most studies emphasized the links between biodiversity loss 
and the risk of pathogen emergence and transmission (Keesing 
et al., 2009; Derne et al., 2011), whereas very few investigated 
the role of parasite diversity, or the size of the 'zoonotic pool', in 
the risk of emergence. 

Investigating parasitic and microbial diversity in bat species 
may help to identify species that are reservoir sources of a greater 
diversity of pathogens and to understand factors influencing this 
richness, particularly those related to life or ecological traits. The 
determinants of parasite diversity in wildlife may be linked to their 
biogeography (e.g., latitude, distribution area), ecology (e.g., den- 
sity, migration), life-history traits (e.g., longevity and fecundity) 
or immunity (e.g., white blood cell counts, spleen size, immune 
gene diversity) (Morand and Poulin, 2000; Wegner et al., 2003; 
Nunn et al., 2003b; Guernier et al., 2004; Ezenwa, 2004; Bordes 
et al., 2007; Simkova et al., 2008; Pedersen and Grieves, 2008; 
Turmelle and Olival, 2009; Bordes and Morand, 2011; Nunn, 
2012; Luis et al., 2013). Some of these determinants can also influ- 
ence the probability of contact of bats with humans and thus the 
potential risk of contamination and transmission; a widely distrib- 
uted bat species living at high density has a greater probability of 
repeated contacts with humans than a species living in low density 
with a restricted distribution. This may be particularly true for syn- 
anthropic species that appear to be generalist in their ecology and 
rich in the parasite diversity they harbour (Herbreteau et al., 201 2). 
Finally, a question rarely investigated about species richness is: "is 
there any correlation between microparasite (viruses, parasitic 
bacteria, protists, fungi) richness and macroparasite (helminths 
and arthropods) richness in bats?" (Bordes et al., 2008; Turmelle 
and Olival, 2009). The diversity of microparasites (e.g., virus and 
bacteria) may then depend on the diversity of macroparasites 
(e.g., helminths) through the activation and maintenance of differ- 
ent pathways of the immune system (Bordes and Morand, 2011; 
Ezenwa and Jolles, 2011). A positive correlation may suggest that 
hosts with high macroparasite diversity also harbour high viral 
and bacterial richness. 

We focused our study on Southeast Asia (SEA), a hotspot of bio- 
diversity and ElDs with pandemic potential (Myers et al., 2000; 
Jones et al., 2008; Coker et al., 201 1 ). It is also a natural laboratory 
to study the evolutionary history (Guillen et al., 1997) and the 
impact of high human environmental pressures (Sodhi and 
Brook, 2006; Clements et al., 2006; Stibig et al., 2007; Wilcove 
et al., 2013). For our study, we defined parasite species richness 
(PSR) as the total number of parasite species such as micropara- 
sites as well as macroparasites identified in a given host (Poulin 
and Morand, 2004; Bordes et al., 2007) at the regional scale of 
SEA. In the present study we aimed at investigating the likely fac- 
tors that may explain PSR in bats from SEA. From the literature we 



compiled information on parasitic and infectious agents in SEA 
bats found, as well as information available on their life-history 
and ecological traits. Then, we tested potential factors that may 
explain the whole pool of parasite diversity in bats using phyloge- 
netic comparative analyses and model selection. 

We tested hypotheses related to (1) influences of the size and 
shape of the geographic distribution (index of the fragmentation 
of bat distribution, defined as the ratio of the surface of the distri- 
bution area to the edge length); (2) the size of bat colonies and the 
number of breeding seasons. We hypothesized that a large area 
should favor the accumulation of parasites and that an increase 
of fragmentation of bat distribution, potentially linked with 
increased habitat diversity, lead to an increase of the overall para- 
site diversity. We also hypothesized that bats species living in large 
colonies with a high number of breeding seasons will support a lar- 
ger pool of parasite species due to large susceptible populations. 

2. Materials and methods 

2.1. Data on bats and their parasites 

Information on bats and their parasites were compiled from the 
literature. Only Southeast Asian countries were selected: Brunei; 
Cambodia; Indonesia; Lao People's Democratic Republic; Malaysia; 
Myanmar; Philippines; Singapore; Thailand; Timor; and Vietnam. 

A total of 292 species bats inventoried in SEA were included in 
the database. Several variables were documented for each host 
species: distribution size; distribution shape (border edge, as a 
measure of the fragmentation of bat distribution area estimated 
by the ratio of the distribution circumference to the distribution 
area size); breeding seasons per year; colony size and average adult 
body mass (Table 1 ). The distribution shape ranges from 0 (cylindri- 
cal and compact distribution area) to 1 (fragmented distribution 
area) (Fortin et al., 2005). Information was gathered from several 
sources: IUCN Red List (http://www.iucnredlist.org/); Wilson and 
Reeder «Mammal Species of the Worlds (http://www.bucknell. 
edu/msw3/); and Harrison Institute (http://www.sc.psu.ac.th/batdb/ 
index.asp). In addition, published articles were obtained through 
match searches using Google Scholar (http://scholar.google.com) 
and the Web of Knowledge (http://apps.weboflmowledge.com/). 

Data on parasite species richness (PSR), which is the number of 
parasite species described infecting a given host species, were 
obtained through searches on Web of Science. To identify relevant 
information binomial names of every Southeast Asian bat species 
were combined with related terms (parasite, ectoparasite, endo- 
parasite, virus, bacteria, fungi, helminths, pathogen and disease). 
All resulting abstracts and available full texts were examined. We 
examined 964 publications from 1959 to July 2012 and identified 
637 species of parasites in SEA. 

In addition we used specific databases for SEA, such as the 
database of Armed Forces Research Institute of Medical Sciences 
(http://www.afpmb.org/content/welcome-literature-retrieval-sys- 
tem); the Liverpool database (http://www.zoonosis.ac.uk/eid2) for 
additional records of viral richness; and the British Natural History 
Museum database (http://www.nhm.ac.uk/research-curation/ 
research/projects/host-parasites/database/index.jsp) for helminths. 

2.2. Sampling effort and investigation effort 

We measured sampling effort, or ascertainment bias, in two 
ways. First, we defined the 'sampling effort' as the sample size or 
number of individual bats tested for a given parasite or microbe 
screened (positive or negative). Second, the 'investigation effort' 
is the number of publications about a given parasite or microbe 
for a given host species, which represents the research effort for 



Table 1 

List of bat species with information on parasite richness (endoparasitic helminths, ectoparasitic arthropods, virus species, virus genus and virus family), sampling effort (SE Para = number of individual bats investigated for a parasite, SE 
virus = number of individual bats investigated for a virus), investigation effort (Pub para = number of publications citing a parasitic of microbial agent in a given bat species. Pub bats = number of publication of a bat species), body mass 
(g), distribution size (km 2 ), distribution shape (area range/circumference), colony size (0 = §10; 1 = §100; 2 = §1000; 3 = >1000), roosting behavior (0 = solitary; 1 = depend; 2 = gregarious), breeding seasons (1 = one season; 2 = two 
seasons; 3 = all year long) and IUCN status (0 = least concern; 1 = near threatened; 2 = vulnerable). 



Bat species 


Endoparasite 


Ectoparasite 


Virus 


Virus 


Virus 


SE 


SE 


Pub 


Pub 


Body 


Distribution 


Distribution 


Colony 


Roost 


Breed 


IUCN 








sp. 


genus 


fam 


Para 


Virus 


para 


bats 


mass 


size 


shape 


size 




season 


status 


Acerodon leucotis 
















1 


2 


400 


12832 


0.1601 


1 


2 




2 


Aethalops alecto 




5 












5 


3 


23.1 


686808 


0.0195 


0 


0 




0 


Aselliscus stoliczkanus 




1 








1 




1 


8 


6.1 


1368853 


0.0069 


2 


2 


1 


0 


Balionycteris maculata 




3 








50 




5 


4 


13.5 


574183 


0.0161 


0 


2 




0 


Cheiromeles torquatus 


1 


10 












15 


6 


181.5 


1338608 


0.0162 


2 


2 




0 


Chironax 




4 












4 


1 


18 


178954 


0.0298 


0 


2 


2 


0 


melanocephalus 


































Coelops frithii 


2 














2 


1 


3.7 


2505073 


0.0159 


1 


2 




0 


Cynopterus brachyotis 


1 


11 


6 


3 


3 


93 


79 


35 


54 


32.8 


2715326 


0.0145 


0 


1 


2 


0 


Cynopterus horsfieldii 


1 


4 


2 


2 


2 




2 


13 


1 


59 


1495518 


0.0153 


2 


2 


2 


0 


Cynopterus sphinx 






1 


1 


1 


81 


81 


1 


128 


47.7 


6489257 


0.0088 


0 


2 


2 


0 


Dobsonia exoleta 




1 












2 


0 


249 


181117 


0.0468 


0 


2 




0 


Dyacopterus spadiceus 




2 












2 


3 


80 


447874 


0.0174 


0 


2 




1 


Emballonura monticola 




5 












5 


0 


35 


1701771 


0.0204 


1 


2 


1 


0 


Eonycteris major 
















1 


8 




737545 


0.0113 


1 


2 






Eonycteris spelaea 


3 


20 


5 


4 


4 


246 


44 


33 


25 


58.4 


3546214 


0.0199 


2 


2 


3 


0 


Eptesicus dimissus 




1 












1 


1 




4087 


0.0673 










Eptesicus serotinus 


2 




1 


1 


1 




1 


27 


128 


25 


12129638 


0.0087 


0 


1 




0 


Glischropus tylopus 


1 


1 












2 


3 


4.5 


2466827 


0.0117 


0 


2 


1 


0 


Harpyionycteris 




1 












1 


1 


120.3 


172958 


0.0722 


0 


1 




0 


whiteheadi 


































Hesperoptenus 


1 














1 


4 


6.4 


463170 


0.0220 


0 


2 




0 


blanfordi 


































Hipposideros armiger 


5 


11 


7 


5 


4 


113 


113 


40 


29 


47 


4496220 


0.0093 


0 


1 


1 


0 


Hipposideros bicolor 




9 








74 




10 


15 


47.7 


1594983 


0.0187 


2 


2 


1 


0 


Hipposideros cervinus 












14 




2 


3 


9.3 


1978116 


0.0241 


1 


2 




0 


Hipposideros 




6 












6 


5 


4.3 


2214839 


0.0124 


0 


2 


1 


0 


cineraceus 


































Hipposideros diadema 


1 


7 


2 


2 


2 




2 


10 


9 


39 


3340330 


0.0261 


3 


2 


1 


0 


Hipposideros galeritus 




4 












6 


2 


7.9 


2358398 


0.0166 


0 


2 


1 


0 


Hipposideros larvatus 




15 


5 


5 


5 


156 


143 


30 


14 


20.4 


3306719 


0.0119 


2 


2 


1 


0 


Hipposideros lylei 




1 












1 


3 


36 


1018679 


0.0116 


1 


2 


1 


0 


Hipposideros pomona 


1 




2 


2 


2 


388 


388 


9 


9 


6.5 


2855362 


0.0092 


2 


2 


1 


0 


Hipposideros pratti 


1 














1 


9 


37.5 


2251718 


0.0062 


1 


2 




0 


la io 






1 


1 


1 




1 


3 


50 


53.2 


1512964 


0.0078 


1 


2 




0 


Kerivoula hardwickii 




1 








15 




2 


7 


3.9 


4070770 


0.0141 


0 


0 




0 


Kerivoula minuta 




1 








4 




2 


2 


2.2 


144409 


0.0282 








1 


Kerivoula pellucida 




1 








2 




1 


5 


6 


1221402 


0.0170 


0 


2 




1 


Macroglossus minimus 




6 


1 


1 


1 


3 


1 


7 


10 


16.5 


3613476 


0.0285 


0 


1 




0 


Macroglossus sobrinus 




2 








19 




2 


5 


22.9 


2449534 


0.0120 


0 


1 


3 


0 


Megaderma lyra 


3 


7 












10 


113 


50 


6070972 


0.0082 


0 


1 


1 


0 


Megaderma spasma 




7 












10 


5 


26.9 


4082464 


0.0193 


0 


1 


1 


0 


Miniopterus magnater 






5 


3 


2 


189 


164 


50 


6 


13.4 


3179606 


0.0128 


2 


2 


1 


0 


Miniopterus medius 




11 












11 


0 


9.5 


1258546 


0.0167 


2 


2 




0 


Miniopterus pusillus 






6 


2 


2 


474 


425 


27 


10 


8.4 


1113321 


0.0166 


2 


2 


1 


0 


Murina aurata 






1 


1 


1 




1 


1 


2 


5.1 


314718 


0.0262 


0 


1 




0 


Murina cyclotis 




1 












1 


7 


7.6 


2032608 


0.0133 


0 


2 




0 


Murina suilla 




1 








4 




1 


4 


3.8 


1451218 


0.0142 


0 


1 




0 


Myotis adversus 




1 












2 


19 


47.7 


244183 


0.0352 


0 


2 




0 


Myotis chinensis 


1 




1 


1 


1 




1 


3 


6 


27.8 


2272978 


0.0059 


0 


0 




0 


Myotis hasseltii 




1 












1 


0 


10 


953379 


0.0173 


0 


1 


1 


0 
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Table 1 (continued) 



Bat species 


Endoparasite 


Ectoparasite 


Virus 


Virus 


Virus 


SE 


SE 


Pub 


Pub 


Body 


Distribution 


Distribution 


Colony 


Roost 


Breed 


IUCN 








sp. 


genus 


fam 


Para 


Virus 


para 


bats 


mass 


size 


shape 


size 




season 


status 


Myotis horsfieldii 




6 


2 


1 


1 


8 


8 


9 


3 


8 


2576798 


0.1454 


0 


1 


1 


0 


Myotis muricola 












1 




1 


5 


3.8 


4668165 


0.0132 


0 


1 


1 


0 


Nyctimene minutus 




1 












1 


0 




26128 


0.0580 








2 


Penthetor lucasi 




10 








10 




13 


2 




1296841 


0.0126 


2 


2 




0 


Philetor brachypterus 






1 


1 


1 


13 


13 


1 


0 


8.2 


1325800 


0.0242 


1 


2 




0 


Pipistrellus abramus 






5 


2 


2 


203 


203 


6 


36 


4.5 


3854526 


0.0133 


1 


2 




0 


Pipistrellus javanicus 




1 












1 


4 


4.2 


3727958 


0.0168 


1 


2 




0 


Pipistrellus stenopterus 




7 












7 


4 


15.5 


644892 


0.0176 


2 


2 




0 


Ptenochirus jagori 






1 


1 


1 


16 


16 


1 


6 


85 


272675 


0.0569 


0 


0 


2 


0 


Pteropus giganteus 


2 




5 


3 


3 


184 


184 


9 


96 


1021 


4021279 


0.0065 


2 


2 




0 


Pteropus hypomelanus 


1 


3 


4 


4 


3 


689 


689 


22 


37 


438 


527562 


0.0625 


1 


2 


1 


0 


Pteropus vampyrus 


2 


3 


2 


1 


1 


475 


475 


34 


52 


911.3 


1953294 


0.0283 


1 


0 




1 


Rhinoiophus 




1 








9 




1 


3 


11.6 


1083857 


0.0199 


1 


2 




0 


acuminatus 


































Rhinoiophus affinis 


6 


8 


2 


2 


2 


292 


266 


25 


13 


12.9 


5534663 


0.0095 


2 


2 


2 


0 


Rhinoiophus lepidus 




6 












6 


9 


7.4 


3526223 


0.0076 


0 


2 


1 


0 


Rhinoiophus macrotis 


2 


3 


4 


2 


2 


51 


46 


9 


12 


5.6 


2348803 


0.0093 


0 


2 




0 


Rhinoiophus 


1 


1 








7 




2 


6 


6.7 


1547550 


0.0095 


2 


2 




0 


malayanus 


































Rhinoiophus marshalli 


1 














3 


5 


6.3 


656136 


0.0099 


0 


2 




0 


Rhinoiophus pearsonii 


5 




3 


3 


2 


95 


95 


8 


2 


13 


3191107 


0.0076 


1 


2 


1 


0 


Rhinoiophus pusillus 


1 




1 


1 


1 


6 


6 


3 


15 


3.5 


3989176 


0.0144 


1 


2 


1 


0 


Rhinoiophus sinicus 


1 




6 


3 


2 


1993 


1993 


28 


10 


11.4 


2177356 


0.0097 


0 


1 




0 


Rhinoiophus stheno 


1 


6 












9 


4 


6.3 


2216461 


0.0129 


2 


2 




0 


Rhinopoma 


4 














5 


13 


31 


7183340 


0.0060 


2 


2 


1 


0 


microphyllum 


































Rousettus 




9 


3 


3 


3 


51 


51 


14 


9 


64.5 


4308473 


0.0217 


2 


2 




0 


amplexicaudatus 


































Rousettus leschenaultii 




4 


4 


3 


3 


16 


16 


8 


7 


79.8 


6798300 


0.0093 


2 


2 




0 


Saccolaimus 




3 












3 


14 


57.5 


1969538 


0.0321 


0 


2 




0 


saccolaimus 


































Scotophilus kuhlii 


7 


10 


7 


6 


6 


284 


284 


44 


12 


21.6 


5254460 


0.0173 


1 


2 


1 


0 


Tadarida jobensis 




1 












1 


0 


25 


2349741 


0.0127 








0 


Tadarida mops 




6 












7 


9 


21.9 


616082 


0.0147 








1 


Tadarida plicata 




11 


4 


4 
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a given parasite investigation (usually corresponds to the number 
of pathogens' screened). Parasite species richness (PSR), viral rich- 
ness, investigation effort and sampling effort were log transformed 
in order to stabilize variance. 

Residual values of the linear regression between log of PSR and 
sampling effort were calculated. We investigated residual varia- 
tions among bat hosts in order to identify bat species that harbour 
more pathogens than expected by the regression model (i.e., bat 
species with positive residual values), or conversely, to identify 
bat species that harbour fewer pathogens than the model predicts 
(negative residual values). This approach can be used to target taxa 
for future pathogen discovery, as was proposed by Herbreteau 
et al. (2012) for rodents and rodent-borne diseases. Of the 292 spe- 
cies of bats inventoried in SEA, information on parasite species 
richness and sampling effort was obtained for 41 bat species, 
whereas information on viral species richness and sampling effort 
was obtained for 33 bat species. 

2.3. Comparative analyses on parasite species richness 

Phylogenetically related species tend to share common charac- 
teristics and are not independent observations (Harvey and Pagel, 
1991). Comparative analysis using the independent contrasts 
method allows studying relationships between species traits by 
limiting false statistical results due to phylogenetically pseudo- 
replications (Type 1 and II errors) (Felsenstein, 1985). We com- 
puted a working phylogenetic tree of bat species investigated in 
this analysis using phylogenetic trees from published studies (see 
Supplementary data). 

Of the 292 species of bats inventoried in SEA, we obtained com- 
plete information on 81, for which we documented all factors 
(independent variables in the subsequent analyses) including a 
fully resolved phylogenetic tree. Analyses were then performed 
on 20 bat species for viruses, 17 bat species for endoparasites 
(helminths) and 28 bat species for ectoparasites. 

We calculated independent contrasts for each of the investi- 
gated variables with the package APE (Paradis et al., 2004) imple- 
mented in R (R Development Core Team, 2008). Independent 
contrasts were calculated for three groups of parasites (viruses, 
endoparasitic helminths and ectoparasite arthropods) for which 
we had a full set of explanatory variables: sampling effort or inves- 
tigation effort; bat body mass; bat distribution range size; distribu- 
tion shape (fragmentation); bat colony size; number of breading 
seasons; and bat gregarious behavior. Parasite species richness 
(PSR), viral richness, investigation effort, sampling effort, range, 
body weight were log transformed in order to stabilize variance. 
Distribution shape was transformed using arcsine of square root 
transformation. To confirm the proper standardisation of contrasts, 
the absolute values of standardised contrasts were regressed 
against their standard deviations (Garland et al., 1992). Then con- 
trasts were analyzed using standard multiple regressions, with all 
intercepts forced through the origin (Garland et al., 1992). 

We selected the model using a backward procedure and due to 
potential co-linearity among variables we performed a Principal 
Component Analysis in order to select variables. We used Akaike's 
Information Criteria (AIC) to select the best models. 

3. Results 

3A. Identification of highly parasitized bat species 

Linear regression between sampling effort and PSR showed a 
positive relationship (P< 0.0001, R 2 = 0.54). Species were ordered 
by residuals of the sampling effort regression (Fig. 1A). Macroglos- 
sus minimus was identified as the host species with the largest 



positive residual value, carrying more parasite species than 
expected by linear regression between species richness and sam- 
pling effort, whereas Rhinolophus acuminatus appeared to host 
fewer parasites than expected (i.e., highest negative residual). Sim- 
ilarly, the regression between virus diversity, a subset of PSR, and 
sampling effort significantly explained 46% of the virus diversity 
variation (P < 0.0001 ). Species were ordered by residual values 
(Fig. IB). Hipposideros armiger carried the greatest number of 
viruses above predictions by the linear regression model, whereas 
Cynopterus sphinx hosted the fewest viruses. 

3.2. Correlation of parasites types per host species 

We found one positive correlation between ectoparasite, endo- 
parasite and virus richness using raw data (Table 2). Ectoparasite 
species richness was weakly positively correlated with virus spe- 
cies richness in bat hosts (P = 0.04). The relationship between ecto- 
and endo-parasites was not significant (P = 0.06); nor was the rela- 
tionship between endoparasite and viral species richness 
(P=0.16). 

3.3. Multivariate analysis 

We performed a Principal Component Analysis on potential 
determinants of parasite and viral species richness of bats (Fig 
2): bat body size; number of publications for each bat species; 
number of publications about parasites of each bat species; bat 
range distribution; distribution shape or fragmentation; colony 
size, breeding seasons; roosting site; and ICUN status. The two first 
dimensions accounted for 43.1% of the variance. IUCN status 
appeared correlated with bat distribution shape (fragmentation 
of distribution area). The number of publications on bat species 
was related to both bat body mass and bat distribution range, 
whereas the number of publications on parasites was related 
to bat distribution range. Bat roosting site was related to colony 
size. 

3.4. Comparative analyses using independent contrasts 

The best model explaining virus richness included distribution 
shape, size of bat colony and investigation effort (Table 3). Both 
distribution shape and colony size were negatively correlated with 
viral richness (Table 4). Bat species with fragmented distributions 
and living in large colonies harbour less virus than bat species with 
continuous distributions and living in small colonies. 

The best model explaining ectoparasite richness identified bat 
body mass, colony size and investigation effort (Table 3). Body 
mass was negatively correlated with ectoparasitic arthropod 
species richness, whereas colony size was positively, but not 
significantly, correlated with ectoparasite species richness 
(Table 4). 

Finally, investigation effort was positively, and the distribution 
shape negatively, correlated with endoparasitic helminth species 
richness (Table 4). 

4. Discussion 

This study is the first investigation of the likely determinants of 
parasite and viral species richness in SE Asian bats, using a large 
data set on parasitic and microbial diversity and ecological traits 
of their bat hosts. The results of comparative analyses showed that, 
in addition to sampling effort, factors related to fragmentation of 
distribution (distribution shape), colony size and body mass seem 
to explain the variability of parasitic (endoparasitic helminths and 
ectoparasitic arthropods) and viral species richness. 
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Hipposideros pomona 
Cynopterus sphinx 
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Rhinolophus acuminatus 




Hipposideros armiger 
Scotophiius kuhlii 
Cynopterus brachyotis 
Eonycteris spelaea 
Minhpterus pusillus 
Rousettus leschenaultii 
Hipposideros iarvatus 
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Myotis horsfieldii 
Tadarida plicata 
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Pteropus hypomelanus 
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Hipposideros pomona 
Pteropus vampyrus 
Cynopterus sphinx 




Fig. 1. Distribution of residual values from the linear relation between (A) Parasite species richness (PSR) and sampling effort and between (B) viral richness and sampling 
effort. 



4.1. Relationship between parasites 

We found a positive and significant correlation between ecto- 
parasite and virus species richness in bat species in our data set. 
This correlation may be related to the vector-borne transmission 
of some viruses by ectoparasites (e.g., Aznar-Lopez et al., 2013). 
For example, ectoparasites, including highly adapted and 



host-specific bat flies (Order Diptera), may play a role as vectors 
for bacteria and protozoa in bats, e.g., Bartonella (Morse et al., 
2012) and possibly Hepatocystis (Garnham, 1951; Olival et al., 
2007). However, the feeding behavior of ectoparasites (blood- 
sucking or not) was not documented as there are few publications, 
but this information could be added in our analyses when 
available. Also due to a lack of data, it was not possible to test 
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Table 2 

Pearson-moment correlations between richness values of ectoparasite, endoparasite 
and virus species in bats (with P-value). 



Endoparasite 
species richness 



Ectoparasite Viral species 

species richness richness 



Endoparasite 

species richness 
Ectoparasite species 

richness 
Viral species 

richness 



0.45 (0.06) 



0.35 (0.16) 
0.47 (0.04) 



cd 

<M 

g 



o 



o 
d 



in 
d 



rang 5 




-1.0 -0.5 0.0 0.5 

Dimension 1 (24.37%) 



1.0 



Fig. 2. Principal Component Analysis performed on the following potential 
determinants of parasite species and viral diversities: bat body size (weight); 
number of publications for each bat species (pub species); number of publications 
of parasites of each bat species (pub para); bat range distribution (range); 
distribution shape or fragmentation (fragmentation); colony size (colony); breeding 
seasons (breeding); roosting site (roosting); and IUCN status (status). 

the correlation between viral richness and bacterial richness. An 
immunologically-driven mechanism of viral and bacterial diversity 
in bats may explain their apparent resistance to virus-induced dis- 
eases and would be worth testing further. For example in humans, 
herpes virus latency confers a surprising resistance to infection 
with bacterial pathogens (Barton et al., 2007). 



4.2. Sampling effort, parasite and viral diversity 

Sampling and investigation effort were positively correlated 
with parasite or viral species richness, as observed in many com- 
parative studies (Walther et al., 1995; Poulin, 1995; Guegan and 
Kennedy, 1996; Nunn et al., 2003a; Turmelle and Olival, 2009; 
Luis et al., 2013). Based on this pattern, we used residual values 
of the linear relationship between PSR and sampling effort as a 
way to draw attention to species with greater than expected num- 
bers of parasites (Herbreteau et al., 2012). 

The association found between sampling investigation (i.e., 
number of publications investigating parasites) and bat distribu- 
tion area, using multivariate analysis, confirms that parasitologists 
mostly screen for parasites of common bat species (i.e., living over 
a large distribution area). 

The ordination of residual values of PSR obtained from the lin- 
ear regression between PSR and sampling effort identified bat spe- 
cies with higher values than expected by the linear regression 
model: M. minimus; Penthetor lucasi; and Rousettus leschenaultii, 
etc. Of the 41 bat species studied, half of species with positive 
residual values are known to carry emerging or potentially emerg- 
ing viruses: Nipah virus; Australian bat lyssavirus; Phnom Penh bat 
virus; Kaeng Khoi virus; and Coronaviruses. It may reveal a high 
level of virus screening in bat species of the same genus. For exam- 
ple, the genus Rhinolophus was considered as a reservoir of a huge 
diversity of bat-SARS-like Coronaviruses in both Asia and Europe 
(Wang et al., 2011; Balboni et al., 2012). Rhinolophus was the target 
of many investigations for Coronavirus discovery (e.g., Lau et al., 
2005; Li et al., 2005; Tang et al., 2006). However, disparities in 
sampling effort among species of a same genus were noted. Rhinol- 
ophus sinicus, Rhinolophus malayanus and R. acuminatus seem over- 
investigated and an increase in sampling effort may not improve 
pathogen discovery, whereas other species of the genus need more 
investigation (e.g., Rhinolophus macrotis, Rhinolophus pearsonii, Rhi- 
nolophus pusillus and Rhinolophus affinis). 

1. Similar ordination of residual values performed for virus species 
richness showed that three bat species have higher residual val- 
ues: H. armiger; Scotophilus kuhlii; Cynopterus brachyotis. Of the 
33 species studied, half of positive residual values belong to 
families Hipposideridae, Vespertilionidae and Pteropodidae. 
Genera belonging to these families are known to cany viruses 
in the families Paramyxovidae and Coronaviridae. For example, 
the genus Pteropus is considered as the natural reservoir of the 
Paramyxovirus, Nipah virus (De Jong et al., 2011; Field and 



Table 3 



Comparison of models used to explain viral richness, ectoparasite and endoparasite species richness in Southeast Asian bats using independent contrasts. The initial model 
included the following variables: sampling effort (number of hosts sampled for a given virus: hsv); investigating effort (number of publications on parasites for a given host: 
pubn); host body mass (body mass); colony size (colony); breeding seasons (breeding); geographical distribution range size (range); and shape of the geographical distribution 
(fragmentation). Models are ranked from the least to the most supported according to corrected Akaike information criteria (AIC). 


Dependent variables 


Model ranks 


AIC 


Viral richness 


Fragmentation + colony + pubn 
Fragmentation + colony + breeding + pubn 
Range + fragmentation + colony + breeding + pubn 
Range + fragmentation + colony + breeding + hsv + pubn 
Body mass + range + fragmentation + colony + breeding + pubn 


31.0 
31.1 
32.6 
34.3 
36.2 


Ectoparasite species richness 


Body mass + colony + pubn 

Body mass + colony + breeding + pubn 

Body mass + fragmentation + colony + breeding + pubn 

Body mass + range + fragmentation + colony + breeding + pubn 


98.4 
99.7 
100.8 
102.3 


Endoparasite species richness 


Fragmentation + pubn 

Body mass + fragmentation + pubn 

Body mass + fragmentation + breeding + pubn 

Body mass + fragmentation + colony + breeding + pubn 

Body mass + range + fragmentation + colony + breeding + pubn 


41.3 
42.7 
44.6 
46.6 
48.6 
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Table 4 

Best model explaining viral richness, ectoparasite and endoparasite species richness in Southeast Asian bats using independent contrasts (initial model with investigating effort, 
host body mass, colony size, number of breeding seasons, geographical distribution size and fragmentation of the geographical distribution) using the A1C criteria (with 
SD = standard deviation of the slope). 



Dependent variables 


Independent variables 


Slope (SD) 


P 


R 2 , F(P) 




Virus richness 


Fragmentation of distribution 
Colony size 

Investigating effort (number of publications on parasites) 


-24.0 (6.8) 
-1.1 (0.2) 
3.6 (0.3) 


0.003 

<0.0001 

<0.0001 


0.87, F 3 , 16 


= 41.77 (<0.0001) 


Ectoparasite richness 


Body mass 
Colony size 

Investigating effort (number of publications on parasites) 


-2.8 (0.9) 
0.5 (0.4) 
8.7 (1.0) 


0.004 

0.18 

<0.0001 


0.76, F 2 , 25 


= 38.89 (<0.0001) 


Endoparasite richness 


Fragmentation of distribution 
Investigating effort 


-72.1 (22.0) 
3.3 (0.7) 


0.006 
<0.001 


0.64, F 2 ,i 4 


= 12.61 (0.0007) 



Epstein, 201 1 ) and was well investigated in SEA. In contrast, it 
seems that several species need more investigation, including 
H. armiger and all the other species with strongly positive resid- 
ual values. This prioritization method highlights the over-inves- 
tigation of some bat species of the family Pteropodidae such as 
C. sphinx, Pteropus hypomelanus, Pteropus vampyrus, Ptenochirus 
jagori, Balionycteris maculata, whereas others need more sam- 
pling efforts for parasites (e.g., P. lucasi, C. brachyotis, Eonycteris 
spelaea). 

Our analysis of residual values between sampling effort and 
parasite and viral richness among bat species can be used to target 
hosts for cost-effective pathogen discovery and also to identify 
hosts that have been well sampled for pathogen discovery. The 
genus Pteropus represents one third of the species with positive 
residual values, i.e., with greater number of parasite and viral spe- 
cies richness than expected by the linear correlation with sampling 
effort. However, species of the genus Pteropus represent 40% of all 
species included in our prioritization analysis and more balanced 
taxonomic sampling in the future will improve the representative- 
ness of analyses similar to ours. Only a few bat species from SEA 
were well documented and thus integrated to the analysis (41 
for PSR and 33 for viruses). As the number of individual hosts 
tested is often missing in published papers, the number of publica- 
tions related to a given parasite in a given host species may be used 
for future residual analyses. This information has the advantage of 
being more easily documented and was correlated with PSR in our 
study. 

4.3. Ecto- and endoparasite species richness 

Ectoparasite species richness was negatively correlated with 
body mass and positively, but not significantly, correlated with col- 
ony size of bat species. A negative correlation of ectoparasite spe- 
cies richness with bat body mass was unsuspected as it was found 
positively correlated in bats from South America (Bordes et al., 
2008). However, using multivariate analysis, body mass appeared 
to be negatively associated with colony size and roosting behavior 
(Fig. 2) and these two last factors could explain ectoparasite spe- 
cies richness. The statistical power of our comparative analysis 
was limited by the number of SEA bat species investigated for ecto- 
parasites. Moreover, our results may suggest that social behavior 
should be defined in a more complex way in order to better under- 
stand factors affecting contact rates, transmission and diversity of 
parasites (Altizer et al., 2003). Consistent with Bordes et al. (2008), 
the distribution area did not seem to explain ectoparasite species 
richness in bats of Southeast Asia. It differs from other studies on 
PSR in terrestrial mammals (Krasnov et al., 2004; Lindenfors 
et al., 2007), where it is assumed that widely distributed mammal 



species accumulate parasites as they increase their chances of con- 
tact with parasite species (Lindenfors et al., 2007). 

Endoparasite richness was negatively correlated with the distri- 
bution shape, an index of the fragmentation of bat distribution. 
One hypothesis is that endoparasites need specific environmental 
conditions to survive and infect bats, mostly linked to the diversity 
of their epidemiological cycles including intermediate hosts and 
external environment stage. Moreover, parasites are known to suf- 
fer from habitat fragmentation that destabilized host-parasite 
interactions (Kruess and Tscharntke, 2000). A bat carrying a wide 
diversity of endoparasites may need a variety of environmental 
conditions to maintain this richness. 

4.4. Viral richness 

The best model explaining viral richness identified investigation 
effort, distribution shape and colony size (both negatively corre- 
lated). The negative correlation observed for colony size does not 
correspond to what would be predicted by classical epidemiologi- 
cal models where a large colony should harbour more viruses 
because of a greater pool of susceptible (larger critical community 
size). One explanation for this finding may be that our measure of 
colony size does not consider the structural complexity of bat com- 
munities. Some bat species live in close proximity with other spe- 
cies. For example, R. leschenaultii can share its roosting site with 
Miniopterus schreibersii, Rhinolophus rouxii, Rhinolophus lepidus 
and other fauna such as rock pigeon (Korad and Gaikwad, 2008). 
An important variable to take into account is then the size of the 
whole community of bat species and the rates of interaction among 
these species. A large colony size does not also equate with high 
density in term of epidemiological transmission. The social com- 
plexity of bats within a colony may limit the overall contacts 
between individuals leading to a decrease of parasite diversity, 
such as observed for the ectoparasite species richness in rodents 
(Bordes et al., 2007). However, information on bat sociality or bat 
community structure is mostly lacking for the investigated bats 
of SEA, and is an area for future research. 

4.5. Effect of the bat distribution shape 

Our analysis highlights the significantly negative effect of bat 
distribution shape on endoparasite and viral species richness. The 
negative correlation between the distribution shape and viral or 
endoparasite species richness corroborate results from Turmelle 
and Olival (2009), who found greater viral richness in bat species 
with more genetically structured populations. However, the rela- 
tionship between genetic diversity and distribution shape (i.e., frag- 
mentation) is not well known but should be a negative as an 
increase in fragmentation should reduce interactions between bats 
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and thus increase genetic divergence. This hypothesis remains to be 
tested as more information on genetic diversity of bats from SEA 
becomes available. Although the distribution shape is not an index 
of the fragmentation of habitat, it may reflect a border effect of bat 
species' distribution area and thus could be considered as a bat spe- 
cies characteristic. This index may originate from the evolutionary 
biogeographical history of each species (expansion and reduction of 
population range) (Hampe and Petit, 2005). An increase of distribu- 
tion shape may reveal species vulnerability, characterized by low 
population size and patchy distribution leading to a decrease in par- 
asite diversity. The importance of the distribution shape draws 
attention on environmental factors that may affect border edge, 
particularly in the face of the ongoing global change. 

Parasite richness could be considered as a zoonotic pool involv- 
ing a risk for human health. According to this view the negative 
effect of bat distribution shape on endoparasite and viral species 
richness could be consider as protective for humans. However, 
fragmentation, due to changes in land use by humans, increase 
rates of contact between humans and animals (Plowright et al., 
2011). This contact may be a critical factor underlying spillover 
(Keesing et al., 2009). Moreover, habitat fragmentation may favour 
emergence as the decline of PSR reduces interspecific competitive 
interactions that could benefits some pathogens. 

Future studies should examine consequences of forest/habitat 
fragmentation impacts on bat distribution and prevalence of ElDs 
(Murray and Daszak, 2013; Olival et al., 2013). Very few studies 
have investigated the effect of distribution shape or fragmentation 
on parasite diversity. Our results are an important first step to 
understand links between PSR and other host and environmental 
traits that may influence disease emergence. Thus, future studies 
should examine effects of habitat fragmentation on the distribu- 
tion shape (with better information on bat distribution) and PSR 
of bats, which may be an important determinant of viral richness 
and PSR more generally. 

4.6. Targeting bats as source of potential viruses 

We are aware that the results of our comparative study were 
strongly dependent on the quality of the published literature, 
and are based purely on correlation analyses. Ordination of resid- 
ual values of the linear relationship between viral richness and 
sampling effort has helped prioritization of several species belong- 
ing to families Hipposideridae, Vespertilionidae and Pteropodidae 
for further viral screening (see above). Ordination of the bat distri- 
bution shape could be used as a way to target the bat species as 
source of potential emerging viruses. Among the first ten species 
that are ordered according to their distribution shape (see Table 1 ) 
appear bat species of the family Pteropodidae (such as C. sphinx, 
Pteropus giganteus or R. leschenaultii), H. armiger (Hipposideridae) 
and some bat species belonging to Vespertilionidae (e.g. Myotis 
chinensis, R. pearsonii, Eptesicus serotinus). However, much work 
should be done to confirm that these species should be targeted 
and to use this kind of comparative analysis as a tool for predicting 
bat species as potential sources of emerging viral diseases. 

5. Conclusion 

We examine determinants of PSR (endoparasite, ectoparasite 
and viral) in SEA bats. First, we used model residuals as a novel pri- 
oritization method to target bat species for a cost-effective patho- 
gen discovery. Secondly, results of comparative analyses suggest 
that the distribution shape is a significant determinant of PSR, as 
well as colony size for viruses and distribution size for ectopara- 
sites. Several potential mechanisms may explain the correlation 
with bat distribution shape through a border effect, although the 



life-cycle and ecology of these parasites may be of importance, 
with opposite trends in different parasite groups. Environmental 
habitat and distribution of bats seem to play a central role in shap- 
ing species diversity of parasites and viruses. Ongoing global envi- 
ronmental change is affecting the distribution of bats through the 
modification of borders, i.e., by decreasing the distribution size and 
increasing the area shape. Our results suggest that changes in bat 
distribution shape will alter parasite diversity, with a decrease 
for endoparasites and viruses. Our study is an important first step 
in understanding parasite and pathogen species richness in bats 
from an emerging disease hotspot. However, accurate investigation 
of emerging risks will require additional information on the role of 
ecological changes in PSR and contact rates between bats, humans 
or their domestic animals. Moreover, our study highlights the defi- 
ciency of ecological and parasitological data on bat species. Of the 
292 species inventoried in SEA only 81 species could be included in 
our comparative analysis. 
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